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Numerical exact diagonalization is the ultimate method of choice in order to discuss static, dy- 
namic, and thermodynamic properties of quantum systems. In this article we consider Heisenberg 
spin-systems and extend the range of applicability of the exact diagonalization method by showing 
how the irreducible tensor operator technique can be combined with an unrestricted use of gen- 
eral point-group symmetries. We also present ideas how to use spin-rotational and point-group 
symmetries in order to obtain approximate spectra. 
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I. INTRODUCTION 

The knowledge of energy spectra of small magnetic 
systems such as magnetic molecules is indispensable for 
the (complete) understanding of their spectroscopic, dy- 
namic, and thermodynamic properties. In this respect 
numerical exact diagonalization of the appropriate quan- 
tum Hamiltonian is the ultimate method of choice. Nev- 
ertheless, such an attempt is very often severely restricted 
due to the huge dimension of the underlying Hilbert 
space. For a magnetic system of N spins of spin quan- 
tum number s the dimension is (2s + 1)^ which grows 
exponentially with N. Group theoretical methods can 
help to ease this numerical problem. A further benefit 
is given by the characterization of the obtained energy 
levels by quantum numbers and classification according 
to irreducible representations. 

Along these lines much effort has been put into the de- 
velopment of an efficient numerical diagonalization tech- 
nique of the Heisenberg model using irreducible tensor 
operators, i.e. employing the full rotational symme- 
try of angular momenta»i^^ A combination of this 
meanwhile well established technique with point-group 
symmetries is not very common since a rearrangement 
of spins due to point-group operations easily leads to 
complicated basis transformations between different cou- 
pling schemes. A possible compromise is to use only 
part of the spin-rotational symmetry (namely rotations 
about the z-axis) together with point-group symmetries 6 - 
or to expand all basis states in terms of simpler prod- 
uct statesi 7 -^ To the best of our knowledge only two 
groups developed a procedure in which the full spin- 
rotational symmetry is combined with point-group sym- 
metries. O. Waldmann combines the full spin-rotational 
symmetry with those point-group symmetries that are 
compatible with the spin coupling scheme, i.e. avoid 
complicated basis transforms between different coupling 



schemes^ Sinitsyn, Bostrem, and Ovchinnikov follow a 
similar route for the square lattice antiferromagnet by 
employing D4 point-group symmetry ! 11 ' 12 This already 
establishes a very powerful numerical method. 

In this article we show how the irreducible tensor oper- 
ator technique can be combined with an unrestricted use 
of general point-group symmetries. The problem, that 
the application of point-group operations leads to states 
belonging to a basis characterized by a different coupling 
scheme whose representation in the original basis is not 
(easily) known, can be solved by means of graph the- 
oretical methods that have been developed in another 
contexti 13 : 14 We discuss how this methods can be im- 
plemented and present results for numerical exact diago- 
nalizations of Heisenberg spin systems of unprecedented 
size. 

Having these methods developed we also discuss ideas 
of approximately obtaining energy spectra of so-called 
bipartite, i.e. non-frustrated, antiferromagnetic spin sys- 
tems. The idea is to perform numerical diagonaliza- 
tions in the orthogonal Hilbert subspaces characterized 
by spin and point-group quantum numbers using only a 
restricted but carefully chosen basis subset. We demon- 
strate how this idea works for archetypical spin systems 
such as bipartite or slightly frustrated spin rings. The ad- 
vantage compared to alternative approximate methods 
such as Density Matrix Renormalization Group 15 ' 16 : 17 
(DMRG), Lanczos, 18 or Quantum Monte Carloi^i 2 ! 
(QMC) techniques is, that one obtains many energy lev- 
els together with their spectroscopic classification which 
can be of great use for the discussion of Electron Param- 
agnetic Resonance (EPR), Nuclear Magnetic Resonance 
(NMR), or Inelastic Neutron Scattering (INS) spectra. In 
this respect our idea can provide a valuable complement 
to the already established approximate methods. 

The article is organized as follows. In Sec. [Til we ex- 
plain the idea of a combined usage of spin-rotational and 
point-group symmetry. Section UTTI provides examples for 
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full diagonalization studies. Our approximate diagonal- 
ization scheme is introduced in Sec. [IVl whereas Sec. [V] 
provides example calculations on bipartite systems. The 
paper closes with a summary. 



II. THEORETICAL METHOD 

A. Irreducible tensor operator approach 

The physics of many magnetic molecules can be well 
understood with the help of the isotropic Heisenberg 
model with nearest-neighbor coupling. The action of an 
external magnetic field is accounted for by an additional 
Zeeman term. The resulting Hamiltonian then looks like 



H = -J2 Ji3*{i) ■ g(i) + 3Msg • B . 



(1) 



The sum reflects the exchange interaction between single 
spins given by spin operators s at sites i and j. For 
the sake of simplicity we assume a common isotropic g— 
tensor. Then the Zeeman term couples the total spin 
operator S = Y^i=i s (*) to the external magnetic field 

B. A negative value of Jy refers to an antiferromagnetic 
coupling. 

For the following discussion an antiferromagnetic 
nearest-neighbor exchange coupling of constant value 
J < is assumed (which can easily be generalized) , then 
the Heisenberg part can be written as 



^Heisenberg = ~J ^ £W ' SO') 



(2) 



<i,3> 



where the summation parameter < i,j > indicates the 
summation running over nearest-neighbor spins counting 
each pair only once. Since the commutation relations 



-^Heisenberg 7 S 







(3) 



hold it is possible to find a common eigenbasis {|^)} of 
-HHeisenberg, § 2 and S z . We denote the corresponding 
eigenvalues as E v , and M v . Due to spin-rotational 
symmetry, Eq. ([3]), the eigenvalues of the Hamiltonian 
([1]) can be evaluated (later) according to 



E U (B) = E v + gfi B BM v , 



(4) 



where the direction of the external field B defines the 
z-axis. 

Calculating the eigenvalues here corresponds to finding 
a matrix representation of the Hamiltonian and diagonal- 
izing it numerically. A very efficient and elegant way of 
finding the matrix elements of Eq. ([2]) is based on the use 
of irreducible tensor operators. Apart from its elegance 
it drastically reduces the dimensionality of the problem 
because it becomes possible to work directly within the 



subspace Tt(S, M = S) of the total Hilbert space TC char- 
acterized by quantum numbers S and M = S; for typical 
dimensions compare for instance Ref. I22I . 

The calculation of matrix elements of the given Hamil- 
tonian using irreducible tensor operators is compulsorily 
related to the application of the Wigner-Eckart-theorem. 
The Wigner-Eckart-theorem 



(aSM\T^\a' S' M') = 
(-l) s - M (aS||T (fe) ||c/ S') 



S 

-M 



k S' 
q M' 



(5) 



states that a matrix element of the g-th component of 
an irreducible tensor operator T( fc ) of rank k is given 

by the reduced matrix element (a 5||T( fe ) | \a' S') and a 
factor containing a Wigner-3J symbol^ 

It should be emphasized that the reduced matrix ele- 
ment is completely independent of any magnetic quan- 
tum number M. The basis in Eq. (|5|) is given following 
the well-known vector-coupling-scheme. The quantum 
number a within the ket \a S M) refers to a set of inter- 
mediate spin quantum numbers resulting from the cou- 
pling of single spins s to the total spin quantum number 
S. In order to apply the Wigner-Eckart-theorem it is 
necessary to express the Heisenberg Hamiltonian in Eq. 
([2|) with the help of irreducible tensor operators. There- 
fore the single-spin vector operators s (i) can be seen as 
irreducible tensor operators of rank k = 1 with compo- 
nents q = — 1, 0, 1. The relation to the components of the 
vector operators is given by 



4( 



s x ±is v 



(0) 



Writing the Heisenberg exchange term as a tensor prod- 
uct of the single-spin irreducible tensor operators results 



-^Heisenberg 

^3 j T^({k l },{k m }\k i = k j 



(J) 



<hj> 



T^ ) is a zero-rank irreducible tensor operator depend- 
ing on the sets {h}, I = 1, . . . , N, which give the ranks 
of single-spin irreducible tensor operators and {k m }, 
m = 1 , . . . , A — 1 , which refers to the ranks of intermedi- 
ate irreducible tensor operators. In a successive coupling 
scheme within a system of N spins an irreducible tensor 
operator of this kind would look like 



T^Uh}, {k m }) = {... {{s^(l) ® s ik ^(2)} ikl) 



.(fc 3 



(3)} 



} 



(A0> 



(0) 



(8) 



The notation T^°\{ki}, {k m }\ki = kj = 1) corresponds 
to the situation in which the ranks of all single-spin tensor 
operators are zero except those at sites i and j which are 
tensor operators of rank 1. 
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The set {k m } results from the chosen coupling scheme, 
for example of the form of Eq. (|8|) , with known ranks of 
single-spin tensor operators taking into account addition 
rules for spin quantum numbers of the vector coupling 
scheme like k\ = \k\ — • • • , k% + k%. 

After writing the Heisenberg Hamiltonian as a sum of 
irreducible tensor operators the matrix elements within 
a basis of the form \a S M) can be calculated by the 
application of the Wigner-Eckart-theorem. The reduced 
matrix elements are determined using the so-called decou- 
pling procedure.-' 4 Since the irreducible tensor operator 
T^ fc ) is given as a tensor product of irreducible tensor op- 
erators with regard to a certain coupling scheme (comp. 
Eq. ||gj|), the reduced matrix element (a S\\T^\\a' S') 
can successively be decomposed into a product of single- 
spin irreducible tensor operators and Wigner-9J symbols. 



B. General point-group symmetries 

The use of irreducible tensor operators for the calcu- 
lation of the matrix elements of the Hamiltonian and 
as a result also of the energy spectrum is essential for 
the treatment of magnetic molecules containing many 
interacting paramagnetic ions. Nevertheless, it is some- 
times necessary to further reduce the dimensionality of 
the problem, either because computational resources are 
limited or a labeling of certain energy levels becomes ad- 
vantageous, e.g. for spectroscopic classification^ 4 ^ Such 
a reduction can be done if the Hamiltonian remains in- 
variant under certain permutations of spin centers. Of- 
ten the spin-permutational symmetry of the Hamiltonian 
coincides with spatial symmetries of the molecule, i.e. 
point-group symmetries, therefore the term point-group 
symmetry is used while one refers to the invariance of the 
Hamiltonian under permutations of spins. 

Using point-group symmetries of the system re- 
sults in a decomposition of the Hamilton ma- 
trix (aSM\H\a'SM) into irreducible representations 

rK)(£?) f a group Q whose elements G(R), i.e. the op- 
erators corresponding to the symmetry operations i?, do 
commute with H . 

The symmetrized basis functions which span the irre- 
ducible representations n are found by the application of 
the projection operator 7?( n ) to the full set of basis vectors 
\otSM) and subsequent orthonormalization. The over- 
complete set of basis states {\a S M T^)} spanning the 
n-th irreducible representation T^(Q) is generated by2£ 

V™\aSM) = (x (n) (#))* \ aSM ) > 

(9) 

where l n is the dimension of the irreducible representa- 
tion r( n ), h denotes the order of Q and x^ n \K) is the 
character of the n-th irreducible representation of the 
symmetry operation R. 



Equation ([9]) contains the main challenge while creat- 
ing symmetrized basis states. The action of the opera- 
tors G(R) on basis states of the form \aS M) has to be 
known. Of course, one could expand \a S M) into a linear 
combination of product states \m\ rri2 ■ ■ ■ tun)- Then the 
action of G(R) leads to a permutation of magnetic quan- 
tum numbers m, within the ket \mx m,2 ■ ■ ■ ttin). But, 
the recombination of the symmetry-transformed product 
states into basis states \aSM) will then be much too 
time consuming for larger systems. 

Following Ref. [l0| the action of G(R) on states \a S M) 
can directly be evaluated without expanding it into prod- 
uct states. Suppose there is a certain coupling scheme a 
in which spin operators s(i) are coupled to yield the total 
spin operator S. Generally the action of operators G(R) 
on states \aS M) leads to a different coupling scheme b. 
Now those states which belong to the coupling scheme b 
have to be reconverted into a linear combination of states 
belonging to a. This is technically a rather involved cal- 
culation, and one would not like to do it by hand and 
for every coupling scheme separately. To the best of our 
knowledge it has never been noted or even used that the 
conversion from any arbitrary (!) coupling scheme b into 
the desired coupling scheme a can be well automatized. 
Suppose there is a state \aSM) a belonging to the cou- 
pling scheme a. The action of a arbitrary group element 
G(R) results in a state \aSM)b belonging to a differ- 
ent coupling scheme b. Then the re-expression takes the 
following form 

G{R)\aSM) a = 

SM) aa (a' SM\aSM) b , (10) 

a' 

where a term like a (a' S M\{3 S M)b is known as gen- 
eral re-coupling coefficient. The calculation of gen- 
eral re-coupling coefficients and the evaluation of Eq. 
(fT0|) can be performed with the help of graph theoret- 
ical methods i 13 ' 14 An implementation of these methods 
within a computer program is a straightforward task (fol- 
low directions given in Refs.EIII]) and one can deal with 
any point-group symmetry. 



III. NUMERICAL EXACT DIAGONALIZATION 

In this section we like to present two applications for 
realistic spin systems that can be treated using irre- 
ducible tensor operator techniques and point-group sym- 
metries, but not otherwise. Both systems - cuboctahe- 
dron and truncated tetrahedron - consist of N = 12 spins 
of spin quantum number s = 3/2 (Hubert space dimen- 
sion 16,777,216). The two spin systems, which are re- 
alized as antiferromagnetic molecules - cuboctahedron^S. 
and truncated tetrahedron^, see Fig. [T] for the struc- 
ture - belong to the class of geometrically frustrated spin 




FIG. 1: Structure of the cuboctahedron (l.h.s.) and the trun- 
cated tetrahedron (r.h.s.).— 
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system s 30 ! 31 ' 32 and are thus hardly accessible by means 
of Quantum Monte Carlo. 
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FIG. 2: (Color online) Complete energy spectrum of the an- 
tiferromagnetic cuboctahedron with s = 3/2. 

Figure [2] shows the energy spectrum of the antiferro- 
magnetic cuboctahedron with s = 3/2. This spectrum 
was obtained using only D2 point-group symmetry which 
is already sufficient in order to obtain sufficiently small 
Hamilton matrices. In addition Fig. [3] demonstrates for 
the subspaces of total spin S = and S = 1 that a rep- 
resentation in the full Oh group can be achieved which 
yields level assignments according to the irreducible rep- 
resentations of this group. 

A complete energy spectrum allows to calculate ther- 
modynamic properties as functions of both temperature 
T and magnetic field B. For the cuboctahedron this was 
already done elsewhere^ 2 - Therefore, we like to discuss 
another frustrated structure, the truncated tetrahedron 
which was synthesized quite recently.— In principle this 
geometry permits two different exchange constants, one 
inside the triangles (Ji) and the other between the trian- 
gles (J2), compare Fig. [TJ A practical symmetry for this 
molecule is for instance C^, whereas the full symmetry 
is Td- Figure [3] displays the complete energy spectrum 
for the case J\ = J2 = J ■ The inset of Fig. [4] magnifies 
the low-energy sector. As in the case of many other frus- 
trated antiferromagnetic systems the spectrum exhibits 
more than one singlet below the first triplet^ 



FIG. 3: (Color online) Low- lying energy spectrum of the an- 
tiferromagnetic cuboctahedron with s = 3/2 in subspaces of 
S = and S = 1. The symbols denote the irreducible repre- 
sentations of the Oh group. 




FIG. 4: (Color online) Complete energy spectrum of the anti- 
ferromagnetic truncated tetrahedron with Ji = J2 = J. The 
inset shows low- lying levels in subspaces with S = 0, 1, 2. 



In Fig. [5] we show the zero-field specific heat (top) as 
well as the zero-field differential magnetic susceptibility 
(bottom). The fine structure of the specific heat, which is 
especially pronounced for s = 3/2, results from the low- 
energy gap structure. The sharp peak is an outcome of 
the gap between the lowest singlet and the group of levels 
around the second singlet and the first two triplets, the 
latter being highly degenerate (both nine-fold including 
Af-degeneracy) . This unusual degeneracy of the lowest 
triplets is also the origin of the quick rise and subsequent 
flat behavior of the susceptibility in the case of s = 3/2. 



IV. APPROXIMATE DIAGONALIZATION 

The previous sections demonstrate that numerical ex- 
act diagonalization in connection with irreducible tensor 
operators is a powerful tool to investigate thermodynam- 
ical properties of large magnetic molecules. Nevertheless, 
sometimes the use of total-spin and point-group symme- 
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We label the eigenvalues E^ and eigenstates \(f>\ ' m 
such a manner that 



Jf>h 




FIG. 5: (Color online) Specific heat c(T,B = 0) (top) and 
differential magnetic susceptibility x(T, B — 0) (bottom) of 
the truncated tetrahedron with Ji — J2 = J. 



tries is not sufficient to obtain small enough matrices. 
For such cases we suggest an approximate diagonaliza- 
tion in this section. The approximation is partially based 
on perturbation theory arguments. First ideas along this 
line were already suggested in Ref . l33l . We will generalize 
and largely extend this idea. 

Let's assume that the spin system is described by a 
Hamiltonian H which acts in the Hilbert space TC. Sup- 
pose there is a zeroth-order Hamiltonian Hq and a de- 
composition according to 

Ho + XH' . (If) 



H 



In the case of non-degenerate eigenstates \4>i } of Hq the 
series expansion 



i&> = \4>D + m > ) + ^\€ > ) + 



Ei = E. 



(0) 



\E. 



(i) 



\ 2 E? ] 



(12) 
(13) 



holds for the exact eigenstates \4>i) and corresponding 
eigenvalues Ei. The index i — 1, . . . , n denotes the states 
of the system. The energy eigenvalues and eigenstates in 
zcroth-order result from a (typically simple or even an- 
alytical) diagonalization of Hq within an arbitrary basis 

of n. 



F (0) ^ p(0) 



Vt = 1, 



, n — 1 



(14) 



holds. Now we do not follow conventional perturbation 
theory as it would lead to a successive introduction of ad- 
ditional terms within the series expansion in Eq. (|12[) . i.e. 
terms with increasing order of A. Instead, we diagonal- 
ize the full Hamiltonian H within a reduced set { 1 4>] ) } , 
i = 1, . . . , ^rcd, of eigenstates of Hq, where n ro( j < n is 
referred to as cut-off parameter. The resulting eigenval- 
ues and eigenstates of this approximation are denoted as 
^approx an( j |^ pp x ), Such an approximate scheme is 
always converging since for n re d = n = dim(7i) all ba- 
sis states are incorporated and the diagonalization corre- 
sponds to an exact treatment of the system, i. e. 



E 



l0 ; 



approx \ 



\4>i) V i . (15) 



It is clear that the speed of convergence depends on the 
choice of Hq. 

The speed of convergence will be different for the vari- 
ous states. Since the approximate diagonalization is per- 
formed with the n re d low-lying states of Hq according 

to (|14[) one expects that the low-lying energy levels con- 
verge quickest against their true values. As in perturba- 
tion theory this assumption relies on the hypothesis that 
energetically higher-lying levels do mix into the desired 
low- lying state with decreasing weight. In perturbation 
theory this expresses itself in the second order corrections 

l<4 0) IWf>| 2 

(16) 



E 



(2) 



E 



E 



(0) 



E 



(0) 



(0) 



which decrease with increasing energetic distance E. 
K (0) . In our approximate diagonalization the diagonal 



(0)\ 



E 



(0) 



\E. 



(i) 



E 



(i) 



H'\ 



(17) 
(18) 



and off-diagonal terms ((j>i\H'\(f>y} of perturbation the- 
ory appear as diagonal and off-diagonal matrix elements 
of the reduced Hamilton matrix. Therefore, the approx- 
imate diagonalization includes zeroth and first order by 
definition and all higher orders partially up to the cut- 
off. The inclusion of eigenstates belonging to degenerate 
eigenvalues of Hq poses no problem in our scheme. One 
should only include all eigenstates of a degenerate eigen- 
value into the approximate diagonalization, otherwise the 
convergence is unnecessarily deteriorated. 



Approximate diagonalization based on the 
rotational-band model 



As a zeroth-order approximation Hq of the isotropic 
Heisenberg Hamiltonian © the rotational-band 
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Hamiltonian2i^i26 



DJ 

"m 



n=l 



(19) 



is chosen^ This choice rests on the observation that in 
bipartite antiferromagnetic spin systems the lowest eigen- 
values within subspaces of total spin S follow the Lande- 
ruleM^- i.e. 



£min(S)- E cx S(S+1) 



(20) 



— in Eq. (UHJ) can be seen as an effective 



The prefactor , x 
exchange constant which couples the sublattice spins S 
to the total spin S of the system. The value of D 



D 



1 



N 1 - 1/JV, 



(21) 



is chosen to match the energy of the ferromagnetic 
state of the system described by an isotropic Heisen- 
berg Hamiltoniani^i N s denotes the number of sublat- 
tices which the classical ground state of the system is 
composed of, Nb represents the number of bonds of the 
system. The eigenstates of -Hrb are analytically given in 
the form \S± ... Sjv s S M), which is an enormous advan- 
tage for the following calculations. The corresponding 
eigenvalues are 



DJ 
'~2N 



5(5 + 1) - S * ( S n + K 



(22) 



The spectrum of the rotational-band Hamiltonian con- 
sists of eigenvalues that form parabolas, so-called 
rotational-bands. In the following a rotational-band is 
defined as a set of eventually energetically degenerate 
eigenstates |5i ... Sn s S M) with fixed values of quan- 
tum numbers S n of the sublattice spins. 




FIG. 6: (Color online) Part of the energy spectra of the 
rotational-band Hamiltonian for a antiferromagnetic spin ring 
N = 8, s = 5/2. Seven super-bands are colored. 



Figure [6] shows the spectrum of the rotational-band 
Hamiltonian for a spin ring of N = 8 spins with s = 5/2. 
The lowest bands refer to a sublattice spin configuration 
of 5i = 52 = 4- 5/2 = 10. The next bands result from a 
deviation of one sublattice spin from its maximum value 
5„. m ax = N/N s ■ s. In such a way the whole spectrum 
can be constructed following Eq. (f2"2"| . The eigenstates 
of the rotational-band Hamiltonian are highly degenerate 
due to the many possibilities of combining single spins 
to the sublattice spins S„ and further on sublattice spins 
S„ to the total spin S. 

Figure [HI also shows that the rotational-band spectrum 
is clustered into super-bands (highlighted by color). A 
super-band contains those rotational-bands for which the 
sum of sublattices spin quantum numbers is the same. 
One clearly sees that within the rotational-band spec- 
trum the low-lying super-bands are well separated. 

Inserting Hub into Eq. (fTTj) yields 



H = H 



Hciscnbcrg 



= H 



RB 



H' 



(23) 



as a starting point for an approximate diagonalization. 
With respect to computational resources and due to the 
fact that the eigenstates of H-rb ar e given in the form 
1 5i ... Sjv s S M) the diagonalization is performed in sub- 
spaces TL (5, M = 5) using the irreducible tensor operator 
technique. In addition, point-group symmetries can be 
used for a further reduction of the dimensionality. How- 
ever, only those point-groups can be applied which do 
not alter the sublattice structure, i.e. do not lead to 
rotational-bands that are not included in the approx- 
imate basis set. Then the symmetry operations on a 
state belonging to a certain rotational-band will always 
produce states which belong to the same band. 



V. BIPARTITE SYSTEMS - SPIN RING 
A. Convergence 

In the following we discuss the properties of the pro- 
posed approximate diagonalization for the example of an 
antiferromagnetic spin ring of N = 8 spins with s = 5/2. 
Figure [7] shows the convergence of the energy levels. In 
order to label the levels the full symmetry group D$ of an 
octagon was used. One clearly sees that the convergence 
within the 5 = subspace is fast and smooth (looking 
almost exponential). 

In subspaces of 5 = 1 and 5 = 2 the convergence 
is also fast, but when only few bands are incorporated 
sharp steps can be observed. This is highlighted by two 
arrows in the bottom graph of Fig. [7] The stepwise con- 
vergence is continued in subspaces with 5 > 2 in a very 
regular way. It can be observed that with increasing en- 
ergy within a certain subspace TC(S,M = 5) the steps 
are slightly washed out. The occurrence of the steps de- 
pends on the rotational-band the states are belonging to. 
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FIG. 7: (Color online) Energy levels of an a antiferromagnetic 
spin ring N — 8 with s = 5/2 as a function of the number 
of occupied rotational-bands used for diagonalization in sub- 
spaces 5 = (top), S = 1 (center) and S — 2 (bottom). The 
arrows in the S = 2 subspace refer to the steps within the 
convergence behavior mentioned in the text. The states are 
labeled according to irreducible representations of D$. 



For example, the energy of the lowest state (i.e. the first 
rotational-band) within Ti(S — 2, M = 2) is decreasing 
if 7 bands are incorporated into the approximate diag- 
onalization while the energies of states belonging to the 
second rotational-band are lowered if 8 bands are incor- 
porated, see also discussion in Sec. IV Bl 

In Fig. [5] the convergence of some low-lying eigen- 
states of this spin ring are presented. The convergence 
behaves in analogy to the convergence of the eigenvalues. 
The stepwise convergence in S — 1 becomes obvious. 
Nevertheless, while using only a fraction of basis states 
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FIG. 8: (Color online) Convergence of the eigenstates of an a 
antiferromagnetic spin ring N = 8 with s = 5/2 as a function 
of the number of occupied rotational-bands used for diagonal- 
ization in subspaces 5 = (top) and 5=1 (bottom). The 
states are labeled according to irreducible representations of 
Da. 



(approximately 30% of the basis states within each sub- 
space) the approximate low-lying eigenstates are practi- 
cally converged against the exact eigenstates. In addi- 
tion, it can be seen that states of higher energy converge 
slower than low-lying states. 

We also investigate the convergence for various single 
spin quantum number s. In Fig. [9] the relative difference 
between the approximate energy values and the exact 
values is displayed for various s in the subspace S = 0. 
The levels which have been chosen belong to the first 
three occupied rotational-bands. One clearly sees that 
the approximate diagonalization converges more rapidly 
the higher the single spin is. This is not surprising since 
the rotational-band model (fill)) , which is based on classi- 
cal assumptions, is itself more accurate the larger s is. 



B. Approximate selection rule 

It turns out that the aforementioned occurrence of 
steps can be understood and even be employed for a fur- 
ther reduction of the size of Hamilton matrices. The 
underlying reason is that the full Hamiltonian connects 
states belonging to different rotational-bands with very 
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FIG. 9: (Color online) Relative difference between approxi- 
mate and exact energy eigenvalues for the lowest states of the 
first three occupied rotational-bands for an antiferromagnetic 
spin ring N = 8 with various s and total spin S — 0. Eq refers 
to the energy of the lowest state of the first rotational-band, 
Ei to the lowest state of the second rotational-band and E-j, 
to the lowest state of the third rotational-band respectively. 



since diagonalization requires an enormous numerical ef- 
fort whereas QMC methods scale much more favorable 
with system size for bipartite systems or even frustrated 
systems above a certain temperature. The role of exact 
or approximate diagonalization then would be to use the 
exchange parameters obtained by QMC for the evalua- 
tion of the enersrv snectrum which then can be used e.e. 
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different strength. After having inspected the reduced 
Hamilton matrices of various bipartite systems we arrive 
at the following empirical selection rule. 

The matrix elements (Si, a S 2 , a S M\H\Si,b S 2 ,b S M) 
of the full Hamiltonian between rotational-band states 
are (several) orders of magnitude bigger than all other 
matrix elements if 

|5i, a -5 2 , o |-|Si |i -5 2 ,6|=0. (24) 

Here S\. a and S 2 . a denote the total spins of sublattices 
one and two in (Si <a £2,0 S M\, respectively. Matrix el- 
ements that are not compatible with this rule can be 
neglected which (after a proper rearrangement) results 
in a new block-diagonal structure of the reduced Hamil- 
ton matrix. These blocks are of smaller size and can be 
diagonalized separately. 



C. Application to {Fei2} 

We now apply the approximate diagonalization to an 
existing molecular spin ring^ 9 . that contains 12 Fe 3+ ions 
with s = 5/2. The system can be modeled by an isotropic 
Heisenberg Hamiltonian with antiferromagnetic nearest- 
neighbor coupling J i 39 ' 40 It was theoretically investigated 
in Ref. 21 with the help of QMC methods; the exchange 
parameter was determined to be J = 35.2 K. 

Our intention is to show that it is advantageous to 
combine a stochastic method such as QMC and an ex- 
act or approximate diagonalization. In such a combi- 
nation the role of QMC would be to determine the ex- 
change parameters from thermodynamical observables as 
done in Ref. [5l|. For large systems this is practically 
impossible using exact or approximate diagonalization 




B (Tesla) 

FIG. 10: (Color online) Approximate spectrum of a spin ring 
with N = 12 spins s = 5/2 calculated using 8 bands, D2 
point-group symmetry and the approximate selection rule in 
Eq. (124[) (top). Corresponding magnetization of the system 
(bottom). The dashed red lines refer to experimental data 
of the first three magnetization steps from Ref. with J = 
35.2 K and k B T/J = 0.01. 

Figure [TU] shows the low-energy part of the approxi- 
mate spectrum of the {Fei2} compound modeled by an 
isotropic Heisenberg Hamiltonian. For the approximate 
calculation of the spectrum full rotational symmetry as 
well as D 2 point-group symmetry are used. The calcu- 
lations are performed using 8 occupied rotational-bands 
in the S = subspace and the corresponding number 
of bands in subspaces with S > 0. Overall 21,570,976 
states have been taken into account, which are only about 
15 % of all basis states (dim(W) = 144,840,476). Addi- 
tionally the approximate selection rule given in Eq. (|24|) 
was used in order to reduce the dimensionality of the 
matrices which have to be diagonalized. Figure [TU] also 



displays the magnetization curve which of course can be 
obtained for a bipartite system by QMC as well. The 
magnetization steps^ can be reproduced using the ap- 
proximate diagonalization. 



D. Next-nearest-neighbor coupling — introducing 
frustration 



In the previous parts we demonstrate that the approx- 
imate diagonalization scheme based on the rotational- 
band Hamiltonian yields good results for bipartite, i.e. 
unfrustrated antifcrromagnetic spin systems. We now 
want to investigate how robust the approximate diagonal- 
ization is against the introduction of frustration. To this 
end we study a spin ring with N = 8 and s = 5/2 with an- 
tiferromagnetic nearest-neighbor coupling J = J nn and 
an additional antifcrromagnetic next-nearest-neighbor- 
coupling J nnn which acts frustrating. In a corresponding 
classical system the Neel state (up- down- up- down- . . . ) 
would no longer be the ground state, instead canting can 
occur. One can qualitatively say that with increasing 
Jnnn/ Jnn also the frustration increases. 

Figure [TT] displays the effect of J nnn in the subspace 
H(S = 0, M = S) for the same system that is discussed 
in Fig. [7] for J nnn = 0. The energy gap between the 
ground state and the first excited state decreases with 
increasing frustration. Moreover, the convergence of the 
ground state as well as of excited states becomes slower. 
With Jnnn/ Jnn = 0.4 the convergence is rather poor and 
the quantum mechanical ground state now belongs to the 
irreducible representation B\ of the symmetry group Dg,. 
This means that the true ground state is not the result 
of an adiabatic continuation (A : — > 1 in Eq. (fTTj) ) from 
the ground state of the rotational-band model, which be- 
longs to A%. We just like to mention for the interested 
reader, that this change of the character of the ground 
state constitutes a so-called Quantum Phase Transition; 
in this case for the antiferromagnetic chain with next- 
nearest-neighbor exchange. 

Summarizing, if frustration is only small the approx- 
imate diagonalization still yields good results. More- 
over, the approximate selection rule is also applica- 
ble which is very helpful in calculating the full spectrum 
of the system. 



VI. SUMMARY 

In this work we have demonstrated that the full spin- 
rotational symmetry can be combined with arbitrary 
point-group symmetries. This enables us to obtain ex- 
actly the complete energy spectrum of Heisenberg spin 
systems for so far unprecedented system sizes. Moreover, 
we have outlined a scheme to approximately diagonalizc 
the Hamilton matrix again using the full spin-rotational 
symmetry and point-group symmetries. This approxi- 
mation works well for bipartite antiferromagnetic spin 
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FIG. 11: (Color online) Energy levels of an a antiferromag- 
netic spin ring TV = 8 with s — 5/2 and additional next- 
nearest-neighbor coupling J„n„ as a function of the number 
of occupied rotational-bands used for diagonalization in sub- 
spaces S = With Jnnn/ Jnn = 0.2 (top), Jnnn / Jnn = 0.3 

(center) and Jnnn/ Jnn = 0.4 (bottom). The states are la- 
beled according to irreducible representations of Dg. 



systems. For frustrated systems the quality reduces with 
increasing frustration. How such a scheme can be re- 
fined for frustrated systems will be the subject of future 
investigations. 
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